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Abstract 

The origin of large magnetic fields in the Universe remains currently unknown. We investigate here a mechanism 
before recombination based on known physics. The source of the vorticity is due to the changes in the photon 
distribution function caused by the fluctuations in the background photons. We show that the magnetic field generated 
in the MHD limit, due to the Coulomb scattering, is of the order ICC 49 G on a coherence scale of 10 kpc. We explicitly 
show that the magnetic fields generated from this process are sustainable and are not erased by resistive diffusion. We 
compare the results with current observations and discuss the implications. Our seed magnetic fields are generated 
on small scales whereas the main mechanisms studied in the literature are on scale bigger than 1 Mpc. However, 
compared to more exotic theories generating seed magnetic fields on similar scales, the strength of our fields are 
generally smaller. 
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1. Introduction 


The origin of galactic and inter-galactic large coherence-scale magnetic fields remains largely unknown 031 - 
Galactic fields of micro-Gauss strength and coherence length between 1 kpc and 10 kpc, and inter-galactic mag netic 
fields of strength 10 7 - 10~ 6 G and coherence length between 10 kpc and 1 Mpc, have been observed [JJ, |2, 3,4], 
iS0i 9]. While the above magnetic field measurements are upper bounds, FERMI measurement of gamma-rays 
emitted by blazars seem to provide lower bound of the order of 10“ 15 G in voids 03.111- 

Several mechanisms have been proposed to explain the origin of the magnetic fields on large scales. These models 
can broadly be categorized into early times and late times. In the case of late time models, fields generated in the proto¬ 
galaxies are spilled to the inter-galactic medium IllMQ, while in the case of an early time models |5. J5|- the 
fields are generated during inflationj Ij§, F7,18, Jjj 20], GUT/EW phase transition [2 j} 3 3 l. radiation domination 
or post-recombination era |24, 3 26]. For example, magnetic fields generated from density perturbation have been 
investigated [27,|28, 29], as well as stochastic background of magnetic fields I 3()i . 3ll 3 . 33[ 34 1 and inhomogeneous 
magnetic fields I35L 36]. Both these — early and late times — categories have advantages and disadvantages. For 
instance, the early time models can not generate required magnetic field strength, while the late time models can not 
generate fields with the required coherence length. All the models provide the seed fields which need to be amplified 
by the dynamo mechanism to maintain the observed galactic and inter-galactic magnetic field strength. 

In this work, we focus on the generation of large-scale seed magnetic fields before recombination. The model we 
consider is based on the work of [26] and does not involve any new physics. The source of vorticity is due to the 
changes in the photon distribution function caused by the fluctuations in the background photons. The model studied 
in [261] suffers from several problems. First, the vorticity was under-estimated: the vorticity is a second-order physical 
quantity and some of the second-order contributions were not taken into account. We correct this by expanding 
the Boltzmann equation up to second order. Second, the electrical conductivity during recombination was taken to 
be arising from Thomson scattering whereas we demonstrate that in the temperature range 20 eV < T < 100 eV, 
the conductivity linked to the Coulomb interaction dominates the magneto-hydrodynamics equation that generates 
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magnetic fields (see Eq. ([8])). However, this effect was corrected in [37]. Third, and crucial point, the current inside the 
plasma was over-estimated. Finally, after recombination, the conductivity of the Universe is very high. The magnetic 
flux can thus be considered as frozen ([38]) and the strength of the magnetic fields generated before recombination 
decrease until the collapse of the first structures. This effect was not taken into account in Ref. [2j5|]. We include this 
effect before estimating a dynamo amplification by proto-galaxies [ 12] and galaxies. 

In Section 0, we present in detail the model of generation of seed magnetic field before recombination. In 
Section ([3]). we compute the main parameters of our model, like the vorticity (13.2b and the conductivity (13.3b . to 
obtain in Section (13.51) a theoretical estimation of the strength and the coherence length of the magnetic field seeds 
generated by our model. In Section <j4]), we compare and contrast our model with other models in the literature. 
Finally, in Section (0, we discuss the importance of our work in the light of the recent measurement. To keep the 
continuity, we have included most of the calculations in Appendices. 


2. The model 


We focus on the pre-recombination era in the temperature range from 100 eV to 20 eV. It is important to note that 
E = 13.6 eV is the ionization energy of the hydrogen, hence below this energy the fraction of free electron x e will 
start decreasing slowly until the beginning of recombination (around 0.25 eV) when x e becomes tiny. For the energy 
range 20 eV< T <100 eV, Universe is well described by a hot and dense plasma at equilibrium composed of protons, 
electrons, photons and dark matter. We neglect ionized helium and any gravitational interaction. We also assume that 
dark matter does not interact with the plasma. Because of the global electrical neutrality of the medium, the number 
density of electrons ( n e ) and the number density of protons (n p ) are equal at 0 th order. Let i e [e, /;), where e stands 
for electrons and p for protons. We expand the number density of particles and their velocity up to second order as in 
[39j], so that 
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We define [5 = n e ln y where n y is the number density of photons. CMB and Big Bang Nucleo-synthesis (BBN) 
constraint the value of/i - 6 x 10 10 0 PHI [ 32 th Although electrons and protons are freely moving in the plasma, 


the plasma itself has to be electrically neutral on the whole. 

Unlike Ref. 112611 . we consider two scattering processes involving the charged particles. We include Thomson 
scattering, which depicts the diffusion of photons on free electrons, and Coulomb scattering of free electrons on 
protons. Let V p , V e - and V y denote the velocities of protons, electrons and photons, respectively. Due to the respective 


inertia of each of these particles, we could expect, like in [26], that V p <c V e <c VC, approximate that the protons do 
not move compared to the electrons, and that the induced current, and the seed magnetic fields /i scc( |, are generated 
by the differential motion of electrons inside the plasma, with the current being thus given by J as -en,,v e . But this 
approach is not correct, because thanks to the tight-coupling approximation and the efficiency of Coulomb scattering 
on the temperature range considered, we will have V p ;$ V e - ^ V y . Because of the inaccurate estimation of the 


current, J and the seed magnetic fields generated were over-estimated in 
in the plasma will be 


The correct expression for the current 


J = e{n p v p -n e v e ') (5) 

« e (n&Vpe + 6n pe v) . (6) 

Finally, as the photons are the most mobile particles in the plasma, the bulk velocity v of the plasma can be 
approximated by v r oc J d 3 p V 7 f 7 , where f y is the distribution function of photons. 
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3. Generation of primordial magnetic field 

From (39] and | Appendix A| we have the following equation describing magnetic fields in the plasma 


d,B s Vx (v e xfi) + -V X 
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where P e is the pressure of the electron fluid, m e the mass of the electron, R = p y /n e m e , 5v ye = v y - v e , T ep (resp. 
T r e) the typical time between two collisions of an electron (resp. photons) on protons (resp. electrons); see l3.3l tor 
more details. Il y is the tensor of anisotropic pressure of the photons. As a consequence, we can distinguish two 
ways of generating the magnetic fields, with different source term. One approach is using curl of the current of the 
plasma, and thus its vorticity, as a source of magnetic field (see for example {26]) in a fluid approximation, i.e. in the 
magneto-hydrodynamics framework by taking into account Coulomb scattering only, so that 


djB seed — V X (v X Rseed) ’ V X ./. 

^Coul 


( 8 ) 


where kq 0 u\ = - 77 —- is the conductivity of the medium due to the Coulomb interaction (cf. Section [33l and |Appendix D| 
and v as v e thanks to the tight-coupling and fluid approximations. The second way to generate magnetic fields is by 
considering the contribution of the pressure of the electron fluid, Thomson scattering and the tensor of anisotropic 
pressure (see for example [42]). 
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In this article, we are going to detail the first process and only refer to the second one to compare the two contributions 
and evaluate if galactic and extra-galactic magnetic fields can be generated by these phenomena. 


3.1. Our model 

Magneto-hydrodynamics (MHD) equation for the evolution of the magnetic field due to resistive diffusion is given 
by 0 


<MW = V X (VxB seed ) + -V X J, ( 10 ) 

K 

where k is the conductivity of the fluid. The first term in the RHS of the above equation is non-zero only after the 
generation of the seed field. Hence, the magnetic field evolves due to the resistive diffusion even in the absence of 
plasma flow. Let us take the curl of the current 

V x J = e(ri [V x 6v pe J + [v«] x 6v pe + 6n pe [v x v] + [vdWp,,] x v). (11) 

To have an idea of the strength of this quantity, we will take into account J ~ e6n ep (V x vj only. As a consequence, 
the magnetic field after a time t is given by 

P 1 e5n P 1 n / \ 

Bseed ~ | dr - -G+ I dr-lGp-QJ, (12) 

Jr, K Jt< K 

where Q = V x v is the vorticity of the plasma, = V x v ; is the vorticity of the fluid j, tj is the time at which the 
perturbation enters the Hubble radius and t is the time at which B see d is generated. We are only going to consider the 
first integral in Eq. (Il2l) . It is important to note that the first term in the RHS of Eq. © is a product of two small 
quantities and hence the above estimate of the field is the leading order contribution. We would like to point out some 
interesting features regarding Eq. ©. First, we have assumed that the diffusion time scales, r = (In/A) 2 , are large. 

In other words, if the diffusion time scales are comparable to the process time scales, then the generation of the seed 
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magnetic fields will not be sustainable and our mechanism will be irrelevant to explain the origin of magnetic fields 
in large scale structures (see [44]). In Appendix H we show that the typical time of diffusion in Eq. ([8]) is larger than 


the age of the Universe. Second, the amplitude of the seed field depends on the amplitude of the vorticity generated 
due to the processes and inversely proportional to the conductivity of the processes involved. Third, the processes 
involved in the generation of the vorticity need to to be the same processes one takes into account in the electrical 
conductivity. Lastly, Eq. ([8]) implies that if the vorticity is conserved, the seed magnetic field is also conserved. 

In the rest of this section, we evaluate different parameters (vorticity, conductivity and coherence length) to obtain 
the primordial magnetic field before recombination. 

3.2. Vorticity 

As discussed above, the cosmological plasma before the recombination is in equilibrium. The velocity of the fluid 
(v) is determined by the Boltzmann equation for the photon distribution function f y : 


(d, + Vy.v}fy(t,£,E,p) = / co , i [/i], 

where the RHS includes the Thomson scattering between y and e~. The fluid velocity v is given by 

^hrf d,m/r - 


(13) 


(14) 


where 14 is the particle velocity. Since, vorticity Q is the curl of the velocity, the second-order terms in temperature 
fluctuations can only lead to non-zero vorticity [26]. To see this, let us expand the photon distribution f y = e~ K ’ {k " T(x ' J]> 
about the average plasma temperature (To) i.e.. 
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where f y 0) = e E ! (k « T u\ Only the term quadratic in 5T/To can lead to non-zero vorticity. For detailed calculation, see 


Appendix B 


For the temperature range [lOOeV, 20 eV], vorticity generated due to Thomson scattering is at least 15 orders of 


magnitude larger than the vorticity due to Coulomb scattering. For more details, see Appendix C Taking into account 
all second-order terms, we obtain: 


O ~ 12 x 10 3 



(16) 


where £ y — 1 is the mean free path of the photons and cr T |, = is the Thomson cross-section and 

A the wave-length of the perturbation. ( y is directly linked to the interaction rate Fn, of the Thomson scattering as 
r Th = l/€ y oc T 3 and from CMB observations, we have the constraint y«3x 10 5 . 

Couple of points are worth noting regarding Eq.[l6] first, by including the second-order terms, the estimate of the 
primordial vorticity has been improved by factor 4 as compared to Ref. [26]. Second, the vorticity (O) is inversely 
proportional to the wavelength of the perturbation A. Physically this is related to the fact that photons tend to diffuse 
from the denser to rarer region and since the photons and electrons are tightly coupled the photons tend to carry along 
electrons leading to decrease in the strength of the magnetic field. In the next subsection, we fix the length scale to 
the Silk damping scale. 


3.3. Conductivity 

One of the crucial steps in our analysis is that the process generating vorticity do not necessarily contribute 
significantly to the generation of the seed magnetic fields. As we have shown in Appendix C Thomson scattering 


contributes to the generation of vorticity, however, we show that the Coulomb scattering contributes to the conversion 
of vorticity to seed magnetic field. 
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The conductivity k of the plasma is given by k - J/E. The ratio of the two conductivities, Kjh and Kcoub is given 

byiSlii 


*Th 
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1 (T 9 x (T M ev)“ 5/2 , 


where ln(A) ~ 10 is the Coulomb logarithm, m e = 0.51 MeV/c 2 the mass of the electron and k H the Boltzmann 
constant (for details, see Appendix D| >. For T < 100 eV, 1/tfcoui > 10/ktii- As a consequence, in the temperature 
range considered, the contribution for the conductivity k in the magneto-hydrodynamic limit in Eq. © to generate 
magnetic fields is 


^Coul — hr, 


2 V2Qr k B T) 312 
inj 2 e 2 ln(A) 


( 18 ) 


For details, see also Appendix E 


3.4. Fluctuations and coherence length 

Unlike the generation of the primordial density (scalar) perturbations, the generation of primordial magnetic field 
crucially rests on the coherence length. If the coherence length is small, then the net magnetic field may effectively 
decrease to zero and the generation mechanism is unsatisfactory. 

To obtain the relation, let us consider scalar perturbations of wave-length A, entering the Hubble radius at time f, 
and temperature 7) » 1 eV. These perturbations produce vorticity in the plasma [cf. Eq. (IT6l) ] which eventually leads 
to the generation of the magnetic field. At time t corresponding to a temperature T (in the range 20 eV < T <100 eV), 
seed magnetic field of strength B see d and coherence length A is generated [cf. Eq. ( 1 1 6[ i I. As mentioned earlier, photons 
tend to diffuse from the denser to rarer region and since the photons and electrons are tightly coupled the photons tend 
to carry along electrons leading to the decrease in the strength of the magnetic field. Hence, we set the wavelength at 
the time of generation of the seed magnetic field to be equal to the diffusion length of the photon. In other words, we 
are looking at wave-lengths that are not affected by the Silk damping 0. Under these assumptions, we get: 


Ctly 
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1/2 
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2 x KT 3 ^ 3 , 
0.85 x 4 3 , 


(19) 

( 20 ) 
( 21 ) 


where Am = Tq/( 1 Mpc) and To is the wavelength of the scalar perturbation today. Physically, Am encodes the current 


coherence length of the magnetic field in Mpc. For detailed computation, see Appendix G 


3.5. Seed magnetic field 

Substituting Eqs. m and CCD in Eq. (fl2T> . we have 


dff seed ~ 1 .1(3 ln(A) / m e c 2 \ 1/2 t6T \ 2 1 3 6n ep 
0 ksTf * \k B T) l t) A 4 n 


( 22 ) 


Then, substituting the corresponding expressions with Eqs. (Il9] >. ( 120b and (|2TT > and integrating from 7) to T (cf. 
Appendix G[ >, we can express the strength of 7? seec j either as a function of the current coherence scale Am of the field 
or as a function of the temperature T at which the seed is estimated 
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References 

Coherence length and strength 

[26] 

A < Mpc, B ~ 10~ 23 G 

[29J 

A > 1 Mpc, B » 10- 23 (d/Mpc) 2 G 

[28] 

T ss 10 Mpc and B « 10' 19 G 

[42] 

lOOMpc > A > lOkpc, B < 10~ 3U G 

[34] 

B < 35 nG 

[27] 

A > 1 Mpc, B* 3 x 10- 29 G 

[31] 

A ss 1 Mpc, B < nG 

[32] 

A ss 1 Mpc, B < nG 

[33] 

A ss 1 Mpc, B < nG 

[51] 

A =3 2 Mpc, B Iec = 5 x 10~ 24 G 

[52] 

A « 10 kpc, B iec = 5 x 10~ 23 G 




Figure 1: Summary of different references presenting different mechanisms to generate seed magnetic fields and their 
coherence length. 


At a coherence length of 10 kpc, our model predicts the seed magnetic field strength to be 10 49 G, which is tiny. 


4. Comparison with earlier results in the literature 

Generation of magnetic field with Eq. (0 with the other source term on similar coherence length 10 kpc has been 
studied in Ref. iR^l . where the evolution of the electron-proton-photon plasma without using the long wave-length 
(MHD) is considered. They obtain magnetic strength of B < 10 30 G from small coherence scale to largest possible 
scale /l e [10 kpc; 100 Mpc]. We can see that in both cases, the seed magnetic fields generated are extremely small. It 
will be very difficult to explain the galactic magnetic field with these processes. 

Over the last two decades, there have been other proposals in the literature to generate seed magnetic fields around 
recombination. For instance, non-linear evolution of primordial fluctuations was investigated in |[29]. In this reference, 
the authors obtained the seed magnetic field strength of B « 10 23 (T/Mpc) 2 G with a coherence-length A > 1 Mpc. 
An other work the same topic ([48]) estimated the amplitude of the magnetic field to be B ~ 10~ 27 G at recombination 
on horizon scale. In Ref. [27], the authors obtained magnetic field strength of B ~ 10 29 (/l/Mpc) 2 G with a coherence- 
length ,1 > 1 Mpc. Another process that has been studied is the generation of magnetic fields based on second order 
cosmological perturbation theory. In Ref. j28ll for example, the authors estimated the strength of the magnetic field 
to be B 10~ 19 (T/Mpc) 2 G for coherence length A ss 10 Mpc. As the reader will be able to notice, the above 
models generate magnetic field over a larger coherence length compared to our model. In [49], it was claimed that 
the generation of seed magnetic fields of strength B 10~ 14 at 10 kpc coherence scale could be explained by second 
order perturbation theory but this allegation was later refuted in ||50] and latter studies in [51 ] estimate seed magnetic 
field of amplitude of 10 24 G and coherence length of 1 Mpc. 

However, even if the mechanism described here seems to be inefficient to produce seed magnetic field strong 
enough to explain the magnetic fields in large cosmological structures, [52] follows an analogous approach with a 
different seed of vorticity arising from condensate of axions. They are able to generate seeds of magnitude 10 24 G 
on coherence length scale 10 kpc. As a consequence, we can expect that a modification of the seed vorticity could 
produce magnetic fields whose amplitude are big enough. 


5. Results and discussion 

We discussed here a model to explain galactic and inter-galactic large coherence-length magnetic fields. We 
considered the generation of seed magnetic fields by the vorticity of the primordial cosmological plasma just before 
recombination, at T <100 eV. We demonstrated that the electrical conductivity in the MHD limit is dominated by the 
Coulomb scattering. 
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We have emphasized that electrical resistivity dissipation time scales [44] are several orders of magnitude larger 
than the age of the Universe (Appendix Hi. This result is compatible with Ref. 112711 . This implies that the magnetic 
fields generated in this process are sustainable and are not erased by resistive diffusion. Our analysis also shows 
that magneto-genesis within the frame-work of MHD is consistent. Since the MHD limit corresponds to large wave¬ 
lengths, the smaller wavelength (large frequency) electromagnetic waves do not propagated. The dominating terms 
are those related to diffusion and amplification [44]. 

We have explicitly shown that the strength of the seed magnetic fields generated before recombination is 10 49 G 
on the current coherence scale of 10 kpc. We emphasize that our mechanism generates seed magnetic fields on small 
scales, whereas the main mechanisms studied in the literature are on scale bigger than 1 Mpc. This scale is directly 
due to the wavelength of primordial perturbation entering the horizon during radiation domination. The advantage 
of small scale magnetic fields is that they are easier to amplify with galactic dynamo, as we discuss it later in this 
Section. Compared to more exotic theories generating seed magnetic fields on similar scales, the strength of our fields 
are generally smaller. However, there are still different options to improve the strength of the seed magnetic fields 
especially with regards to the vorticity which is needed to generate them. 

It is difficult to compare our model with current observations of galactic and extra-galactic magnetic fields on kpc 
coherence scale. Based on FERMI observations of Te V sources at z ~ 0.1, lower limit on inter-galactic magnetic fields 
fimin were established to be /( nlm = 10 18 - 10 1 s G at large scales in voids (see Refs. [Il0il53l] ). The observational 
constraints of large-scale magnetic fields in voids give an estimation of the seed magnetic fields. It is important to point 
that some authors have critized the latter tests. In Ref. [54] it was argued that the lack of an inverse Compton GeV 
bump cannot be used as a constraint on the inter-galactic magnetic fields because of plasma instabilities. However, 
the presence of plasma instabilities in TeV blazars beams involve an upper constraint on the inter-galactic magnetic 
fields: B < 10 12 G. Our results are compatible with this constraint. 

In Appendix I we have considered a more realistic process for galactic dynamo amplification. After the first 
stage of galaxy formation our model predicts B ~ 10 47 G. To explain the current observed magnetic fields in the 
galaxies and clusters of the order of 10 6 G coherent on scales of about 10 kpc, our model requires a huge galactic 
amplification. For massive galaxies, like the Milky Way, it is possible to have large amplification, however for small 
and/or lighter galaxies, like dwarf galaxies (sstl , the current dynamo models do not provide large amplification. Also, 
the dynamo mechanism is not efficient in the inter-galactic medium. We have taken into account a dynamo effect as 
realistic as possible, but there is still room for improvement as the different dynamo effects still present numerous 
unknown features. 
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Appendix A. Plasma equations 

From [39], we have the following equations 


d,v e + (v e .V)v e + Hv e = 


VP, e 


-(E + v e x B) 


e~n e n p p y \ 

-Op - V e ) + - Oy - Ve ~ TC-fly) 

*fcoul T ey 4 


- V 0 


(A.l) 


^ -> Vr„ e -> . -> 

d,v p + (\’p.W)vp + Hv p = -- + — (E + VpXB) 


Pp 


e~n e n p m e \ p y 1, n , 

- (v e ~ v p ) + — — (Vy - v p - -Vp.Uy) 

*Coul \m„ T ey 4 


(A.2) 


where {v,. Pi,pi, n,} are the velocity, pressure, energy density and particle density of the fluid i, n y the tensor of 
anisotropic pressure of the photons, 4> the gravitational potential, Kcoui the conductivity of the medium (due to Coulomb 
scattering), r ye the typical time between two shocks of a photon on electrons (Thomson scattering), and E and B the 
electric and magnetic fields in the plasma. Subtracting Eq. (1A. 1 b to Eq. (1A.2I >. we can deduce 


d,5v pe + [(v.v) 5v pe + (dv pc .v)dv] + H5v pe ~ e + ~ j ^ 




m e V 1 Py 


m p J nipiip x e y 


*Coui \m p m, 


Vy ~ v p - -vp.n y 


1 Py 
Wl e Tle Tey 


Ve - Vy + -i4.n y 


Pp Pe 


(A.3) 

The term in j (~ 10' 6 ) for protons can be considered as negligible in front of the analog term for electrons. 

We also neglect the quantity 1 /m p compared to 1 /m e , and v p /m p compared to v e /m e . In addition to the previous 
simplifications, for i e {e,p} 


VP/ = ± 

Pi mi 


vr + t- 


. Vn, 


(A.4) 


At 0 th order, n e ~ n„ so that — * —. With m e /m p ~ 10 3 , then we have H-^H || —[[. As a consequence, 

t n e ftp r Pp p p 


dtdVpe + [(v.V) dv pe + (6\?pe$) 6v] + H6v pe ~ + —(E + v e xB)- - v , pe 

LV ' V ' J Pe m e ' ’ m e Kc out 
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(A.5) 


9 
























Knowing that * 6v pe , 

d, 


f 



/ 


J 

\ 

6v 


f 

+ 

(v.v) 

+ 

.V 

+ H 

[en e 


V / 

[en e 


[en e 

J 



[en e 


VP, 


VTa C / -> , -A 

- + — £ + v,xB ■ 

p e v ' 


1 p r 


2 

e n 


^p^Coul 

1 


en e 


v e - Vy + -Ve.Ti-y 


(A. 6 ) 
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Let us take the curl of the previous equation. Knowing the equation of Maxwell-Faraday Vx£ = -d,B, we thus have 


V x 3, 
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5v 
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+ V x 

(v.V) 
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e ti 
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+ V x 
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v e - v y + -vt.n. 


The terms on the left side of the equation are very small compared to \\-^-d t P\\ because 

-^d,B -^d,B 

m 1 m 1 


kd, 


— » 1 


We can finally write the following equation ruling the behavior of magnetic fields in the plasma 


d,B »Vx (v e xfi)+-V: 


VP„ 


!l e 

v e / 


- Vx 


* p/? c 


Rm e I 1 

V X I —•V^.ily 

eTy e y 4 


where R - p y /n e m e . 


(A.7) 


(A. 8 ) 

(A.9) 


(A. 10) 


Appendix B. Computation of vorticity: upto second order in distribution function 

In this Appendix, we provide detailed computation of the derivation of the vorticity in Eq. (IT6l) by expanding the 
Boltzmann equation for the photon distribution f y up to the second order. The Boltzmann equation for the photon 
distribution f y is given by 


( d t + V.V - Hp.djj + f.dpjf y (t,x,E,p) = / co n[/i], (B.l) 

where V* is the photon velocity. It is important to note that the product of an odd number of 14 is equal to zero whereas 
V k V, = 5 . 

Let us define the operator 'K - d, + L.V - // p.d r ; + P.dp. In our analysis we have ignored the effects of gravity 
(cosmological expansion) as we are working on time scales such that the expansion does not drastically affects the 
results. We also assume that there is no other external force acting on the plasma. As a consequence, 7C = d t + V'.V. 
The velocity of the fluid is given by 


Vk = 



(B.2) 
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For the temperature range considered, the Thomson scattering rate dominates the collision term in the RHS of the 
Boltzmann equation. 


r Th K t 3 oc i y x 

(B.3) 

5/r T h = 3r T hy- 

(B.4) 

3,r T h = 3r T hy- 

(B.5) 


We assume that the elastic electron-photon rate is high and that the integrals over time are dominated by small values 
of time r. We can thus expand Fxh(f — t 2 , x - vt 2 ) as following 


exp 


T Th (t - t 2 , x - vt 2 ) 



dT 2 r Th (t - T 2 , X - VT 2 ) 


r Th (f, x)-T 2 [3,r Th (f, x) + V,<9 ( T T h(f, X)] 


g-Tlf-rhO,:?) x 


1 + ^ T i r Th(f, X) 


dj + d£ 

T ' T 


Appendix B.l. 0 th order: f y x 

At 0 th order, f y 0) = e~ E/T and there is no collision term / co ii[/i] in the Boltzmann equation 


•Kf™ = 0 



As a consequence, the 0 th order, no vorticity is generated in the fluid. 

Appendix B.2. 1 st order: f y » f y 0> + f y X) 

At first order, the collision term is no longer equal zero in the Boltzmann equation 


(B.6) 


(B.7) 


(B.8) 

(B.9) 


Solving Eq. (IB. 1 1L we obtain 
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We expand f y 0) (t - t\,x- vY|) in the same way as we expandedF(f - t\,x- vti) 
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Let us now expand the exponential function in rL and integrate over time 
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Finally, we can compute the fluid velocity at first order v® 
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Appendix B.3. 2 nd order: f y « /* 0) + F 11 + 


y ~ ./y wy wy 

Proceeding in the same way as above, the second order expansion of the Boltzmann equation is given by: 


(7C + r Th )/f = -TC/W 


/3 
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rV 

Jo 


driexp 


- f dr 2 r T h(^ - T 2 ,x- vt 2 ) 
Jo 




Expand F'V? — n, j?— VTi ) as we did earlier for /y U, (f - ri, x - vt\ ) and F(f -T\,x - vt\) 
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Neglecting contributions to the veloctity arising from the terms oc ^ with n > 2, we get. 
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We can deduce from the former equation the contribution to fluid velocity at second order in 1 /T 
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Finally we obtain the following expression for v t2) 
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Appendix B.4. Final expression 
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Due to the symmetry properties, we can notice that the term d,TjT does not contribute to vorticity. Using the fact 
that ly oc 1 !n e « J- 3 , we have 


O, 


tijkdp’k 


8jT 


-24 x r y e ljk ^— x 
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We have the following equation of diffusion in the plasma (see |26H) 



Switching in Fourier space, we get 


lyC r. 

d,T = -yk 2 T 

As a consequence, we can estimate the vorticity 
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By taking into account second order term in fluid velocity, the vorticity we have obtained is 4 times larger than the 
one that was previously obtained in Ref. [26,]. It is important to note that when the temperature decreases the fraction 
of free electrons also decrease, and just before recombination the plasma is no longer at equilibrium and Eq. (HI is 
no longer valid. 


Appendix C. Ratio of Vorticity generated due to Thomson and Coulomb scattering 

The vorticity generated in the Plasma by taking in to account the Thomson scattering is given by: 


Oxh = [---] x dr r x exp 


[...] x e r 


J dr r x exp - J' dr 2 r T h (t - r 2 , x- Vr 2 ) 


whereas the vorticity generated in the Plasma by taking in to account the Coulomb scattering is given by 


^Coui = [-.] x dr r x exp 

Jo 

* [-]x4 2 . 

We can now compute the ratio of the two vorticities 


dr rx exp - J' dr 2 r Cou i (f -r 2 ,x- Vr 2 ) 


(C.l) 
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As a consequence, for the temperature range considered, Oxh » flc; ou i. and Thomson scattering is the dominating 
counterpart for the generation of vorticity in the plasma. 


Appendix D. Ratio of Thomson and Coulomb conductivities 

The equation of motion of an electron in an electric field is the following 


m e 


dVe 

dr 


-eE. 


(D.l) 
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The conductivity k — J/E, where J = -en e V e , of the plasma is given by 

e 2 n e At 


(D.2) 


where At is the typical time between two collisions. The typical time between two electron-photon shocks is (cf. [26]) 

Afxh = (D.3) 

Vj 

where 


4 = c 


3 m e 1 
kgT CTj^Hy 


is the mean free path of the electron and Vj = x/3kfjT/m e is the electron thermal velocity. The typical time between 
two electron-proton shocks is [56] 


A?Coul — tcp 

6 ( nk B T) v2 

n p e 4 ln(A) 

As a consequence, we can deduce the conductivity due to Thomson scattering 
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and the conductivity due to Coulomb scattering a 
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evaluate the ratio of the two conduc- 


^Coul 


= 6 V2< 
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0 mj 1 e 2 n p ln(A) 


To estimate the most important contribution into the MHD equation, let us 
tivities 
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(TMe\) 5/2 

=> *Th = 

Stoul <=> T 

= 250 eV 

(D.ll) 


At T\ = 650 eV, l/Ar T h = 10 /kcou1 and at Tj = 100 eV, l//ecoui = lO/zerh- 


Appendix E. Conductivity due to Coulomb scattering in the MHD equation 

Appendix E.l. Coulomb vs Thomson 

In this Appendix, we show that the Coulomb scattering contributes significantly in the generation of magnetic 
field in the temperature range of our interest. Let us consider the following MHD equations for the two scattering 
processes: 

d,B Th = V x (Lx B-ih) + —V x / (E.l) 

d,S c oul = v x (v x zlcoul) + v x / (E.2) 
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The total contribution from the two processes can be approximately written as 


• (^Th + f^Coul) - V X [V X (z?xh + ^Coul)] + (-1-j V X 7 


*Th *Coul / 


we can define an effective magnetic field 


$ seed - ^Th + ^Coul 


(E.3) 


and an effective conductivity 

11 1 

K Kxh *Coul 


(E.4) 


At T\ - 650 eV, 1 /k-\ \ } = 10/kcouI and at 74 = 100 eV, 1 /kcou1 = 10/x-rh- As a consequence, the contribution for the 
conductivity k in the magneto-hydrodynamic limit in Eq. © is 



KCoui = 6e 0 2 #if^ for T < 100 eV 

u m e ' e l ln(A) 

for T > 650 eV 


. ^coui^Th/ (^Coui + *Th) for T e [100 eV; 650 eV] 
As a consequence, the equation for the generation of the seed magnetic field is 


(E.5) 


djB see d — V X (v X flseed) V X J 

K 

To be more rigorous, an additional term //B see d should be added on the right side of the equation. However, as 
the typical time scale of the generation of the magnetic field is smaller than the expansion time, it will not have much 
impact on the estimation of the field and we can consider it as constant multiplying B sca i 10- 

Appendix F. Expression of the seed magnetic field 


B seed ~ exp(-fcvf) 


fdr^ 

1 Ak 


exp(&:i'T) + constant 


(F.l) 


where kvt ~ 0.015 T/l eV is a small quantity, so we can consider that exp(-kvt) ~ 1 and exp(-£vf) « 1. As a 
consequence, we have 


/Tec;! 


J' dr- (ti5Q. ep + £ldn ep ^ 
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(F.2) 

(F.3) 

(F.4) 


with T(td = Tj, T{t\) = T\ and Titi) = T%. For the temperature range considered, T < 100 eV, we have the following 
rough estimate of the seed magnetic field generated by Coulomb conductivity 


df? seed 


1.7 x (k B T)' 


ffln(A) ( m e c 2 V l2 /dT\ 2 £ 3 y8n e p 
l k eT ) \T) A 4 n 
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Appendix G. Strength and coherence length of B see d from primordial fluctuations 

In this Appendix, for a given physical quantity M, M, refers to the current value of the quantity at Hubble radius, 
M the value at the time of the generation of the magnetic field, and Mq its current value. 


Appendix G.l. Primordial fluctuations 


Let us consider a fluctuation entering the Hubble radius at radiation domination far from 1 eV. The Hubble param¬ 
eter is 158, 26[] 


H7 


27 kpc 



(G.l) 


For T > 20 eV, \/{ y = cr T hti e , hence, 


( y j ~ 30 pc 



(G.2) 


In the radiation dominated era, a oc f l/2 . As a consequence, // 1 = 2t, and we deduce from Eqs. (1G.1I) and (IG.2b 

that 


Cyi , leV 

— ^ 2.2 x 10” 3 x-. (G.3) 

ctf T t 

By definition, the wavelength of the fluctuation entering the Hubble radius is given by T, = cH ~ 1 . As we are 
interested in To, which gives us the wave-length of the fluctuation today and thus the correlation-length of the magnetic 
field today, let us express the temperature T, as a function of T w = To/1 Mpc. We know that the temperature evolves 
as T oc T _1 , so 


li = - 

T T, 


(G.4) 


and y where T y 0> = 2.7255 K is the temperature of the CMB today (see [3^] for the latest evaluation). We 

Ty I 

can deduce that 


Ti = 110 x Aj eV. 


(G.5) 


Appendix G.2. Evolution of the fluctuations 

For the temperature range considered, the seed magnetic fields are generated during radiation dominated era. Now, 
let us focus at the time in which the seed magnetic fields are generated. As a consequence, as the temperature evolves 
as T oc or 1 , 



Using the fact that l y oc T 3 , we have 



(G.6) 


(G.7) 
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In order to keep our fluctuation unaffected by Silk damping, and to have a vorticity strong enough to generate our seed 
magnetic field, we take A ~ £ d where £ d is the diffusion length of the photon. This diffusion length is equal to 158. 2bh 


« = 


— CtCy 

3 r 


(G.8) 


From Eqs. (IG.6b . (IG.7I) and (IG. 8b . we have 


from Eqs. (1G.4I) and ( 1G.9I >. and by remembering that A, = H 1 = 2cf,-, we get 


(G.9) 


'■y.i 


Tj) 6 ctj 


(G.10) 


Appendix G.3. Useful expressions 

From Eqs. (lG.3b . ( IG.5I) and (IG. 1 Ob . we can deduce that 


and from Eqs. (lG.3b and (lG.5b . that 


— 

= 0.015 x A 

r, 

T 

= 1.65 T m 2/3 , 

1 eV 


ty.i 

2 x 10 _5 /1 m . 

Ctf 


1/3 
M ’ 


fG. 11) 


(G.12) 


We finally get three last formula that are used in Section (13.5b to evaluate the strength and the coherence-length of 
the magnetic field seed. From Eq. (IG. 1 Ik we get 


UV V ,2 = 0. 8 5X,« 


T ) 


from Eq. (lG.8b . we have 


Ctly 3 

~ 2 


and using Eqs. (fcTTTb . dGJ2li and (fcTTTh . we obtain 


j) -2xl0' 3 4 3 
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Appendix G.4. Resulting seed fields 

For the temperature range considered, T < 100 eV, we have the following rough estimate of the seed magnetic 
field 


dfiseed 

deseed 

1 G 


Sn e , 


We can estimate — 21 with 


1.7 x(k B T)- 
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fG. 13) 
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5n ep » ^^ V.py \6v yb - ^-Vh.Tiy 
Let us neglect the anisotropic pressure term II y , we are reduced to 


(G.15) 


With the help of (390, we can express v y b. 


6n ep » [p r (5v yi ]. 


(G.16) 


As a consequence, we get the final expression. 


nm p c 2 6 y 
4ahc A 2 


(G.17) 




ep 


136 
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fG. 18) 


Then, substituting Eqs. (1G.4I >. (IG. 6 b and (1G.18I) in Eq. (1G. 13b . integrating between T , and T and knowing that 
Y ^ 1 , we finally have 


^seed 

Tg" 


3.7 x 1(T 52 x T 512 


Appendix H. Typical time scale of diffusion 

Let us now estimate the typical time of diffusion r in Eq. (O to check if the process we are studying is sustainable. 


T (COUl) = poKCoul X (^:) 

6nsj2^ lk B T\ l,2 tk B T\ 1 / A \ 2 
aln(A) \m e c 2 ) \ he ) c\2n) 
= 4.3 x 10~ 4 x x A 2 

= 5.8 x 10 34 x — 

r r'/^ 

1 


(H.l) 

(H.2) 

(H.3) 

(H.4) 


As a consequence, T^l° | l jJ ( , ) eV = 5.8 x 10 27 s. As the age of the Universe is approximately 4.36 x 10 17 s, we do not 
need to worry about the dissipation of our seed magnetic fields for T <12 keV. 
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Appendix I. Amplification of the seed fields 


First, let us consider a seed magnetic field B see d generated in the temperature range considered earlier just before 
recombination. The conductivity of the primordial plasma after recombination is very high, therefore the flux of the 
magnetic field is conserved and the magnetic field evolves as B oc a 2 [38]. This property is taken into consideration 
until z ~ 10. This redshift corresponds to the lowest redshift at which the collapse of the early structures, such as the 
first proto-galaxies, is expected to happen |59[ 60[ 61,62], As a consequence, the magnetic field decreases by 2 orders 
of magnitude until z ~ 10. 

Then, when the first structures collapse, the seed magnetic field gains 4 orders of magnitude using a pre-amplification 
due to adiabatic compression in the pre-galactic medium before galactic formation [12]. The combination of the two 
effects described above in this Section amplifies the seed magnetic field generated before recombination by 2 orders 
of magnitude. 

Finally, for z e [0,10], large-scale magnetic fields are maintained in galaxies and clusters by a dynamo amplifica¬ 
tion. For spiral galaxies, the most popular galactic dynamo model is the aa >-dynamo whereas other theories are used 
for elliptical galaxies and clusters (see 111 [63] for a detailed and critical discussion on the different dynamos). The 
dynamo amplification JK, from an initial field strength B\ at time t \, corresponding to the end of galaxy formation, to 
a galactic field strength Bq at time fo, i.e. our current epoch, is given by 


= — = e r( '° _tl) , 

B i 


(1.1) 


where T is the amplification rate. F is highly dependent on the cosmological model chosen. However, it is usually 
found in the literature that 0.2 Gyrs < T 1 < 0.8 Gyrs (2]. Small values of T 1 , such as T 1 = 0.2 Gyrs, are related to 
processes on small scales, i.e. at intra-galactic scales, and involve fast amplifications. Large values of T _1 , for example 
T 1 = 0.8 Gyrs, are related to processes on large scales, at inter-galactic scales, and trigger slow amplifications. The 
galactic dynamo amplification can thus lie in the interval [10 7 ,10 28 ]. As a consequence, magnetic fields on large 
scales are going to be less amplified than those on small scales. This is coherent with the data, which reveal that there 
not much difference of strength between this two types of magnetic fields. This also guarantee that the magnetic fields 
in the inter-galactic medium are not to strong compared to data. 
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